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We study the dynamics of the East model, comprising a chain of uncoupled spins in a downward- 
pointing field. Glassy effects arise at low temperatures T from the kinetic constraint that spins 
can only flip if their left neighbour is up. We give details of our previous solution of the non- 
equilibrium coarsening dynamics after a quench to low T (Phys. Rev. Lett. 83:3238, 1999), including 



the anomalous coarsening of down-spin domains with typical size d 



t' 



and the pronounced 



'fragile glass'-divergence of equilibration times as = exp(l/r'^ ln2). We also link the model to 
the paste-all coarsening model, defining a family of interpolating models that all have the same 
scaling distribution of domain sizes. We then proceed to the problem of equilibrium dynamics at 
low T. Based on a scaling hypothesis for the relation between timescales and lengthscales, we 
propose a model for the dynamics of 'superdomains' which are bounded by up-spins that are frozen 
on long timescales. From this we deduce that the equilibrium spin correlation and persistence 
functions should exhibit identical scaling behaviour for low T, decaying as g{t). The scaling variable 
is i = (i/f*)"^'"^, giving strongly stretched behaviour for low T. The scaling function g{-) decays 
faster than exponential, however, and in the limit T ^ at fixed t reaches zero at a finite value of 
i. 
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X 



I. INTRODUCTION 

The phenomenology of glassy systems — see e.g. P, 0, 
U 3 for exceUent reviews — has inspired many theoreti- 
cal descriptions and explanations. Experimentally, long 
relaxation times are observed; when these become much 
longer than the observation timescale a glass transition 
is said to occur. Other signatures of glassy dynamics are 
correlation functions that can be fitted by a stretched 
exponential decay law and ageing phenomena where, 
since the system is out of thermal equilibrium, it keeps 
evolving as time goes by and time-translation invariance 
is broken. 

From a modelling perspective the same phenomenolgy 
arises when one studies simple model systems by com- 
puter simulation. Again, relaxation times can outstrip 
the time available to run a simulation and one never ex- 
plores the equilibrium state. 

The long relaxation times in glasses typically show 
a pronounced divergence as the temperature T is low- 
ered and are often fitted experimentally by the Vogel- 
Tammann-Fulcher (VTF) law 



T = ToexphA/(r-To)] 



(1) 



The relaxation time t may characterise, for example, the 
time for a density fiuctuation or an externally imposed 
stress to relax. Although some heuristic justifications 
have been offered Q, for practical purposes VTF is just 
a fit with three parameters tq,A,Tq. For Tq = it re- 
duces to an Arrhenius law. A system for which Tq is 
small, so that one has something close to Arrhenius be- 
haviour, is referred to as a 'strong glass', whereas a sys- 



tem exhibiting large deviations from Arrhenius behaviour 
is called a 'fragile glass'. Generally, Tq is much lower than 
the experimental temperatures so that the mathematical 
singularity in the fit Q is not physically relevant in an 
experiment. From the theoretical point of view, however, 
there has been a long debate over whether Tq might repre- 
sent a true thermodynamic transition temperature which 
would in principle be measurable in the limit of infinitely 
slow cooling. 

Although is popular, it is not the only possibility for 
a fit. For example, the exponential inverse temperature 
squared (EITS) form 



T-roexp(B/T2) 



(2) 



has been proposed as an alternative. This form does 
not exhibit a singularity at any finite T. Experimentally, 
or in a computer simulation, it is difficult to distinguish 
between VTF and EITS behaviour due to obvious lim- 
itations on the longest accessible timescales; both can 
represent the experimentally observed t(T) in many ma- 
terials 7] . Theoretical work is thus essential for clarifying 
whether VTF or EITS might be more appropriate. 

Stretched exponential decay of a relaxation function, 
let us say an autocorrelation g(t), is expressed by the 
Kohlrausch- Williams- Watt law 



g(t)~exp[-(t/r)'' 



(3) 



where the stretching exponent 6 < 1. An heuristic expla- 
nation for this law is that there is a broad distribution 
il(T) of relaxation modes with decay constants r. 



J dr fl{T) exp(- 



-t/r). 



(4) 
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For example if one assumes r2(T) ~ exp(— or) then for 
large t the dominant modes have r = (t/a)^^^ which leads 
to © with b = 1/2. This however leaves the physical 
mechanisms by which such a relaxation time distribution 
would arise unclear. 

One idea proposed to generate a broad distribution of 
relaxation times from an explicit dynamical model was 
that of a hierarchy of degrees of freedom ^ . The different 
levels in the hierarchy then relax in series, the degrees of 
freedom in one level having to wait for the degrees of 
freedom in the level below to reach some configuration 
before they are free to evolve. This latter condition is a 
realisation of a kinetic constraint. 

A more concrete realisation of a kinetic constraint is 
the /-spin facilitated kinetic Ising model introduced by 
Fredrickson and Andersen Q . Here no hierarchical struc- 
ture needs to be posited by hand to generate slow dynam- 
ics. All microscopic degrees of freedom are of the same 
kind, namely non-interacting Ising spins on a simple cu- 
bic lattice in a downwards pointing field. However a spin 
can only flip if at least / nearest neighbour spins are 
pointing up (against the field). For / > 2 this gives rise 
to slow cooperative relaxation, as reviewed in e.g. |lC| . 
The physical motivation for the kinetic constraint in this 
model becomes clear if one thinks of the spins as coarse- 
grained density variables in a supercooled liquid, with up 
and down-spins corresponding to low- and high-density 
regions. The constraint then states that a change in lo- 
cal density is only possible if the overall density in the 
surrounding regions is low enough for the particles to be 
able to rearrange. 

An interesting modification of the Fredrickson- 
Andersen model arises when the kinetic constraints are 
made anisotropic. In the one-dimensional version consid- 
ered here, called the East Model > ^ spin can only flip 
if its left (west) neighbour is pointing up, so that informa- 
tion propagates only to the east. The rates for flipping 
mobile spins are 1 for down-flips and e = exp(— l/T) 
for up-flips, giving a small equilibrium concentration 
c = e/(l + e) of up-spins at low T (see Sec. ^ for de- 
tails). The number of spins which are mobile because 
their left neighbour is up will then also be small, and 
this makes it plausible that the dynamics will slow down 
dramatically as T decreases. 

In this paper we study the glassy dynamics of the East 
model, both in equilibrium and out of equilibrium. We 
focus on low temperatures T. This is the regime that is 
most interesting since glassy features will be most pro- 
nounced; it is also the regime where theoretical studies 
such as ours are most needed, since the extremely long re- 
laxation timescales make numerical simulations difficult 
or impossible. 

Before outlining the structure of this paper, we give a 
brief review of existing work on the East model. Much 
research to date has concerned the equilibrium dynam- 
ics, as encoded e.g. in relaxation functions such as spin 
autocorrelation C(t) or persistence 7^i(i) (see Sec. IIV A1 
for deflnitions). Already when the model was first pro- 



posed [ill it was argued that relaxation timescales should 
remain finite for any T > 0, thus excluding a transition 
to a nonergodic state where relaxation functions fail to 
decay to zero. This has recently been proved rigorously: 
the longest relaxation time, defined as the inverse of the 
smallest decay rate that one would find by full diago- 
nalization of the master equation, is bounded between 
exp(l/2r^ In 2) and exp(l/r2ln2) in the limit of small 
temperatures . The form of the bounds demonstrates 
that the East model has EITS behaviour; we will find 
below that it is the upper bound = exp(l/r^ln2) 
that gives the dominant low-T behaviour of the relax- 
ation times. 

A number of approaches have been used to predict the 
actual shape of the equilibrium relaxation functions. An 
'effective medium approximation' for C{t) leads 

to a self-consistency equation typical of mode-coupling 
approximations (MCA) and predicts a spurious noner- 
godic transition at c = 0.5 (e = 1); effectively the same 
result was later derived using diagrammatic methods |3 ■ 
The version of MCA derived by Kawasaki (l^l also gives 
a spurious transition, at c = 0.2 (e = 1/4). Both ap- 
proximations can therefore only be reasonable at suffi- 
ciently large e, or for short times at smaller e; a com- 
parison with numerical simulations [l^ shows that the 
effective medium approximation is generally more accu- 
rate in these regimes. Improved diagrammatic resumma- 
tions ^3 avoid the prediction of a spurious transition, 
and are quantitatively more satisfactory over a larger 
range of t and e. However, for small e they still pre- 
dict a decay of C{t) that is too fast and too similar 
to an exponential compared with numerical simulations. 
Similar comments apply to the results of 17], where an 
adiabatic approximation together with the hierarchical 
nature of the relaxation processes — discussed in detail 
below — was used to study the analogue of C (t) for finite 
chains. The non-exponential behaviour of C{t) had been 
noticed early on [ll|, but is well fitted by a stretched 
exponential only over a limited time range. We conjec- 
tured in [l^l) and will find below, that the correct scaling 
variable for C{t) at low T is indeed of a 'stretched' form, 
i — (t/t*)"^'"^, but that the relevant scaling function de- 
cays more quickly than an exponential. The stretching 
exponent Tin 2 was also found to govern the decay of the 
up-spin persistence Vi{t) [l9ll20| . 

Recently, interest in the East model has shifted to 
more complicated features of the equilibrium dynam- 
ics, e.g. the existence of dynamical heterogeneities |2l| . 
and to the out-of-equilibrium behaviour. In (2^, non- 
linear relaxation processes after large changes in T were 
simulated. Otherwise, interest has centred on the be- 
haviour after a quench from high to low temperature. 
For a quench to exactly T = 0, the dynamics is exactly 
solvable [23. i2^ . The solution is essentially equivalent 
to that of the corresponding model with isotropic con- 
straints, i.e. the one-dimensional Fredrickson- Andersen 
model with f ~ 1 and in fact also applies to 

a whole family of models that interpolate between the 
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isotropic and anisotropic limits [19|, |20j. For a quench 
to nonzero T, the autocorrelation function C{t,t^) be- 
tween spins at times and t after the quench and the 
corresponding response were simulated in "23]. The cor- 
relation C{t,t^) exhibits plateaux, which we will sec in 
Sec. lIII (!l can be rationalized from the hierarchical nature 
of the dynamics. Non-equilibrium steady states caused 
by applying an external 'drive' have recently also been 
studied, using a 'tapping dynamics' inspired by ideas 
from granular media |27| as well as 'rheological driving' 
designed to model the effect of a shear flow |2^. Fi- 
nally, we mention an interesting two-dimensional spin 
model, which can be mapped onto a system of non- 
interacting defects with kinetic constraints and turns out 
to have low-T behaviour very similar to that of the East 
model 

The paper is organised as follows. After defining 
the East model fully in Sec. |nl we turn to the out-of- 
equilibrium dynamics in Sec. IIIII giving details of the 
results announced in After a quench to low tem- 

perature, the equilibrium concentration of up-spins at 
the new T is small compared to its initial value. Thus 
up-spins are eliminated, essentially irreversibly, and the 
dynamics can be viewed as a coarsening process whereby 
down-spin domains coalesce as up-spins disappear. We 
find the scaling of the rates at which domains of length 
d disappear, which follow a hierarchical pattern so that 
the coarsening dynamics splits into well-separated stages. 
This allows us to find an exact solution for T ^ 0. The 
overall scaling of the rates, as T{d) ~ d^^/"^''^^, also im- 
plies that the system exhibits anomalous coarsening, with 
typical domain sizes increasing with time as d ~ i^'"^. 
Extrapolating to equilibrium domain lengths deq ~ 1/e 
(see below) then gives the dominant divergence of the 
equilibration time as ^ exp(l/r^ In 2). Interestingly, 
in the scaling limit of the out-of-equilibrium dynamics 
(1 <C d ^ f^cq), we find that the domain size distribution 
is identical to that of the paste-all model 32]. We ratio- 
nalize this observation, and provide in Sec. IIII Dl a new 
family of models with identical scaling distributions that 
interpolate between the paste-all and East model limits. 

In Sec. IIVI we shall then explore the equilibrium dy- 
namics. Based on a plausible conjecture regarding the 
scaling of relaxation times with distance at low T, we 
introduce a simplified picture which we refer to as the 
superdomain model. This should capture the essence of 
the equilibrium dynamics and become exact in the T — > 
limit. The model enables us to probe the long timescales 
of the equilibrium dynamics numerically, and leads us to 
a new scaling form for the equilibrium spin autocorre- 
lation function, Clt) = git), with i = (i/i*)^^"^^ This 
shows strong stretching for T ^ 0, but it will turn out 
that the scaling function g decays more quickly than an 
exponential, actually decreasing to zero at a finite value 
of t. Finally, a brief summary of our results and outlook 
to future work is given in Sec. 



II. MODEL DEFINITION 

The model comprises L Ising spins = 0, 1 on a one- 
dimensional lattice with periodic boundary conditions 
(site i = L -I- 1 is identified with site i = 1). The dy- 
namics are defined by the following spin-flip rates 

11^10 with rate 1 

(5) 

10^11 with rate e = exp(-l/r) . 

Thus a spin can only flip if its left neig hbour is pointing 
up (note that in the original paper the mirror image 
of the above definition was used, so that the right neigh- 
bour had to point up for a spin to be able to flip). By 
a rate, say cc, we mean that in a small time dt the event 
happens with probability xdt. It is easy to check that 
the dynamics obey detailed balance with respect to the 
energy function E = equilibrium distri- 

bution corresponds to free spins in a downwards pointing 
field: 

^eq(W)^|exp(^-l^..^=^i|:^. (6) 

It follows that the equilibrium concentration of up-spins, 
c = (si), is given by 



To show that this is the unique steady state we require 
that the (finite) system is ergodic, i.e. that any configu- 
ration can be reached from any other. This is true for all 
configurations except for the one with all spins down; this 
is a configuration which can be neither entered nor left 
and which we therefore ignore. To see that the rest of the 
configuration space is ergodic, note that from any config- 
uration with a non-zero number of up-spins one can fiip 
spins up until the 'all up' configuration is reached. Then, 
to obtain any desired configuration one fiips spins down 
in an ordered way to create the appropriate regions of 
down-spins. 

The basic objects that we use for the description of the 
system are domains. As shown by the vertical lines in 

...0|1000|1|1|10|100|1|1|10|10..., 

a domain consists of an up-spin and all the down-spins 
that separate it from the nearest up-spin to the right. 
(This convention is opposite to that of |3i but leads 
to equivalent conclusions and will be more convenient in 
our treatment of equilibrium dynamics in Sec. lIVI below.'l 
The length d of a domain then gives the distance between 
the up-spin at its left edge and the next up-spin to the 
right. Note that adjacent up-spins are counted as sepa- 
rate domains of length d = 1. In equilibrium, it follows 
from © that the domain lengths are geometrically dis- 
tributed 

Peq(d) = e/(l + (8) 
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with mean 

rfcq = ^ ■ (9) 

As explained, our aim will be to pursue analytical calcu- 
lations where possible, but we explain briefly how sim- 
ulations were carried out. We used a continuous-time 
BKL algorithm, where the time intervals between 
spin flips are sampled directly; a standard Monte Carlo 
algorithm where spin flips are first proposed and then 
accepted or rejected would be much slower at low T . In 
the BKL algorithm, once the time interval to the next 
flip is determined one decides probabilistically which of 
the mobile spins to flip, choosing each with a probability 
Pi proportional to its flip rate (1 or e). A simple method 
for doing this would be to sample a uniform random vari- 
able r on [0, 1] and then go through the mobile spins un- 
til the spin i is found for which X]j=i-Pj <r < Yl!^j=iPi- 
This search for i takes 0{L) steps, however, and in fact 
for large L quickly becomes the most computationally 
intensive part of the algorithm. Instead, we define par- 
tial sums qj^a over all blocks of spins of length 2°" , where 
a = 0, . . . , Hf the total chain length is L = 2'. The qj^a 
can be defined recursively by qj^a=a = Pj ^md 

Qj.a = q2j-l.a-l + q2j,a-l (j = 1, . . . , L/2°) 

(In fact, we work with analogous partial sums for the 
number of mobile up- and down-spins separately, from 
which the qj a can easily be retrieved; this allow us to 
use fast integer arithmetic.) The qj^a can be thought of 
as arranged on a binary tree, and finding the spin i to 
flip becomes a simple walk from the top level a = ^ to the 
bottom level a = 0, branching left or right on each level 
to keep r within the gj^a-values bounding the remaining 
subpart of the tree. This takes ^ ~ InL steps, and the 
qj^a can be updated in similar time once a spin has been 
flipped, giving a reduction in computing time ~ (In L)/L. 
Compared to earlier simulations, e.g. |22| . rather longer 
timescales (up to t = 10^^) can be accessed with this 
method. 



III. NON-EQUILIBRIUM DYNAMICS 
A. Coarsening dynamics 

In this section we consider the dynamics of the East 
model after a quench from equilibrium at some high ini- 
tial temperature T 3> 1 to low temperature T <C 1 
(e 0) At the new temperature the equilibrium 

concentration of up-spins is much smaller than before, so 
that the number of up-spins must decrease in time. Cor- 
respondingly, the typical domain sizes must grow: the 
system coarsens. 

To understand the nature of this coarsening process, 
recall first that the equilibrium concentration of up-spins 
is, from O, c = l/deq = e+0{e^). Hence the equilibrium 



probability of finding an up-spin within a chain segment 
oi finite length d is 0{de) and tends to zero for e — + 0. In 
the limit 



e ^ at fixed d , (10) 



the flipping down of up-spins therefore becomes irre- 
versible to leading order 34]. In terms of domains, this 
means that the coarsening dynamics of the system is one 
of coalescence of domains: an up-spin that flips down 
merges two neighbouring domains into one large domain. 
Note that the limit (|10|l does not extend to include equi- 
librium domain lengths that scale as d ~ 1/e. Therefore 
the irreversible coarsening is a purely non-equilibrium 
phenomenon and does not pertain to the equilibrium dy- 
namics; this is why at low T the former is in fact easier 
to analyse than the latter, as we shall see. 

Irreversible coarsening processes such as the one above 
have been studied in a variety of contexts. In particular, 
if the rate of elimination depends solely on the domain 
size and not on the sizes of neighbouring domains there 
is a very convenient property: during such a process, no 
correlations between the lengths of neighbouring domains 
can build up if there are none in the initial state. The 
proof is an easy generahsation of that of 35] . We include 
it for completeness here. Consider first an initial arrange- 
ment of N domains on a ring. There are A^! ways for this 
system to coarsen, i.e. N possibilities for the first do- 
main to disappear, A^ — 1 possibilities for the second etc. 
Each possibility is weighted with a probability which is a 
function of the rates of disappearance of the A^ domains. 
Summing over the (A^ ~ 1)! possible initial arrangements 
of the A^ domains on the ring implies N\{N ~ 1)1 possible 
'histories' for the dynamics. 

Now note that if the domains are uncorrelated then 
the dynamics becomes equivalent to a mean-field model 
(often referred to as a bag model or an independent in- 
tervals theory): when a domain is eliminated one picks 
a domain at random to be its right neighbour and coa- 
lesces the two. In this dynamics there are at each step 
N{N — 1) possibilities: the factor N comes from which 
domain is eliminated and the factor A^ — 1 from which 
domain is picked to be its neighbour. Therefore for A 
initial domains there are Nl{N — 1)1 possible histories, 
each weighted by the probability for the domains to be 
eliminated in the specified sequence. These possible his- 
tories are in one-to-one correspondence with the possible 
histories for the full model when summed over all possi- 
ble initial arrangements. One also easily verifies that the 
probabilities with which these histories occur, as well as 
the distributions for the times that separate successive 
events within each history, are the same in both cases. 
This proves that throughout the irreversible coarsening 
process the domains are uncorrelated. 
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FIG. 1; Elimination of a domain of size d. Shown are paths 
through spin configurations that traverse the minimum energy 
barrier. The height of the barrier is h{d) and the excited spins 
are underlined in the highest energy configuration(s) along the 
path. 



B. Energy barriers 

We now estimate the typical timescale T~^{d) for the 
disappearance of domains of length d through coalescence 
with their right neighbours. This then defines a typical 
rate T(d) for the elimination of domains. Because domain 
coalescence corresponds to the flipping down of up-spins, 
r{d) can also be defined as follows. Consider an open 
spin chain of length d, with a 'clamped' up-spin (sq = 1) 
added on the left. Starting from the state (sg, si, ■ • • , Sd) 
— 10... 01, T^^{d) is the typical time needed to reach 
the 'empty' state 10 ... 00 where spin Sd has relaxed, i.e. 
has flipped down. Any instance of this relaxation process 
can be thought of as a path connecting the two states. 
Let us call the maximum number of 'excited' spins (up- 
spins except So) encountered along a path its height h. 
One might think that the relaxation of spin Sd needs to 
proceed via the state 11. . . 1, giving a path of height d. 
In fact, the minimal path height h(d) is much lower and 
given by 



h{d) = 



1 for 2"-^ < < 2'^ 



(11) 



where n = 0, 1, . . . 

To get a feeling for the result (|ll|l consider in Fig. 
some small domain sizes. The flgure illustrates that to 
generate an up-spin adjacent to the one to be relaxed 
(which is the first spin of the next domain on the right) 
one can proceed via a sequence of stepping-stones, e.g. for 
d = A one first generates an isolated up-spin in the middle 
of the domain then uses this stepping-stone to generate 
the subsequent excited spins in a similar manner to the 
relaxation of a d = 2 domain. Iteration of this argument 
straightforwardly proves (|11|1 for d = 2" . 

In order to prove Hill) generally we introduce the quan- 
tity l{h, k) which is defined for the open chain discussed 
above. It is the length of the longest configuration that 
contains exactly k up-spins, always counted ignoring the 
fixed up-spin sq, and that can be reached from an ini- 
tially empty chain (again, except for sq) along some 



path of height h. The length of a configuration is de- 
fined here as the position of the furthest up-spin to the 
right. Note that l{h,k) > l{h,k — 1); thus l{h,h) gives 
the longest configuration that can be reached along any 
path of height h. 

Now consider how a configuration realising l{h,k) 
could be constructed. First one generates an isolated 
up-spin as far to the right as possible. This distance is 
l{h, 1). Then one starts from the up-spin just generated 
to generate a second up-spin as far to the right of the first 
as possible. Since the first up-spin remains up, there is 
one less unit of energy to play with to generate the sec- 
ond up-spin and the total distance from the origin will be 
l{h,l) + l{h — l,l). Continuing in this fashion to generate 
all k up-spins, one arrives at 



l{h,k) ^^l{h 



1,1) 



(12) 



To close the set of equations (|12|l we need an expres- 
sion for l(h,l). This can be obtained by invoking the 
reversibility of paths: since the dynamics obeys detailed 
balance, the existence of a path of height h from an initial 
to a final state implies that the reversed path connects 
the final to the initial state and is also of height h. Thus 
l{h, 1) gives the size of the longest domain that can be 
relaxed along some path of height h. It is easy to see that 
the final step in this process requires an up-spin adjacent 
to the up-spin which is to be relaxed. Since, with this 
latter up-spin, we have already used up one unit of en- 
ergy or height, the maximum distance from the origin at 
which such an up-spin could be generated is l{h — l, h — 1). 
Therefore we conclude 



l{h,l) = l{h- l,/i- 1) + 1 



(13) 



with boundary condition ^0,0) — 0. Inserting H13() 
into (|12|1 gives 



l{h, k) — l{h — m,h — m) + k 



(14) 



m=l 



Setting k = h, has solution l{h,h) = 2^-1. Substi- 
tuting this back into (|14|) yields 



h nh—k 



l{h,k) = 2" - 2 



(15) 



In particular the longest domain that can be relaxed via 
a path of height h = n + 1 is given by l{n -1-1,1) = 
l{n,n) -1-1 = 2" from which we deduce (|11() . Related 
results on the number of configurations reachable at or 
below height h can be found in Ref. [s^ . 



C. Solution of the dynamics 

From the result (|ll|l it follows that the coarsening 
dynamics naturally divides into stages distinguished by 
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n = h{d) — 1 = 0, 1, . . . During stage n, the domains with 
lengths 2"^^ < d < 2" disappear; we caU these the 'ac- 
tive' domains. This process takes place on a timescale 
of 0{r~^{d)) — 0(6^"); because the timescales for dif- 
ferent stages differ by factors of 1/e, we can treat them 
separately in the limit e ^ 0. Thus during stage n ac- 
tive domains are eliminated. The distribution of inactive 
domains {d > 2") changes because elimination of an ac- 
tive domain implies coalescence with a neighbouring do- 
main on the right; since the smallest active domains have 
length d = 2"^^ + 1, any new domain will have length 
> 2(2"-i + 1) > 2" and thus be inactive. 

Let N{d,t) be the number of domains of length d at 
time t, N{t) — J2d ^{d, t) the total number of domains 
and P{d,t) = N{d,t)/N{t) the domain size distribution. 
Then for a general process of coarsening by coalescence 
one has, using that there are no spatial correlations be- 
tween domains 

§iN{d,t) = -r{d)N{d,t) - N{d,t)Y,^{d')P{d',t) 



^N{d~d\t)rid')P{d',t) 



(16) 



The first term accounts for disappearance of a domain by 
coalescence with its right neighbour, the second for co- 
alescence with a domain on the left of length d' , and 
the third for creation of larger domains during coa- 
lesence. Summing (|16ll over d one finds dN{t)/dt = 
-Y,^T[d)N{d,t), and thence 



-r{d)P{d, t)+Y^ P{d - d', t)T{d')P{d', t) 

d' 

(17) 

We now apply this general result to the scenario at 
hand. We change to a rescaled time variable appropriate 
to the n-th stage of the dynamics, r = te". To be in 
the n-th stage, t has to obey the restrictions e~("^i) <^ 
t < £-("+!), giving e < T <C e"^ in the limit e ^ 
we can thus let r range over [0, oo]. The rescaled rates 
r(d) = e~"r(c?) are non-zero only for the domains with 
d < 2"; on the other hand, P{d,T) is nonzero only for 
d > 2"^^ since smaller domains have disappeared during 
earlier stages. Thus the sum over d' the r.h.s. of H17|) 
is restricted to 2"^^ < d' < 2", i.e. the active domains. 
Furthermore, if d is an active domain then d — d' < 2"^^ 
and hence the factor P{d — d', t) vanishes for every term 
in the sum, reflecting the fact that active domains cannot 
be created. Eq. lfT7|) thus becomes for 2"'^ < d < 2" 



d_ 



P(d,r) = -f(d)P(d,T) 



(18) 



For d > 2", on the other hand, the first term in (|17|l does 
not contribute and one gets 



J2 P{d-d',t)f{d')P{d',t) (19) 



Combining the last two results one finds that the rates 
drop out, giving for d > 2" 



-i<d'<2" 



d\ 



(20) 

There is one subtlety in this derivation which we have 
so far ignored. Even in terms of the rescaled time r, i.e. 
taking into account only those relaxation processes which 
cross the minimum energy barrier of n, the 'survival prob- 
ability' for a domain not to have coalesced with its right 
neighbour is not a single exponential which can be char- 
acterized by a rate constant T{d). This is because several 
relaxation process with different rescaled rates of 0(1) 
may exist. Thus in 1)181191) one should replace T{d) by the 
negative derivative of the survival probability, resulting 
in an effective, time-dependent rate T{d, t — t') measured 
from the time t' when the domain was created. As ex- 
plained above, however, domains that are active cannot 
be created in the current stage; thus r' must correspond 
to a previous stage, which in the limit e — > implies 
t' 0. Thus H18I19|) are correct if f (d) is replaced by 
r(d, r), and H20() follows as before. Note that, indepen- 
dently of the details of its derivation, equation H2U|I has a 
simple intuitive meaning. It expresses the fact that active 
domains eliminated during the n-th stage coalesce with 
neighbouring domains of random size, forming stable do- 
mains of size d > 2". Equation (|2()|l then determines the 
time dependence of the size distribution of the inactive 
domains for a given time dependence of the distribution 
of the active domains. 

To solve the equation of motion (|20l) we define the 
generating function 



G{z,t)= Pid,r)z' 

d=2"-i + l 



(21) 



and its analogue for the active domains, 

2" 

H{z,t)= P{d,T)z''. (22) 

(i=2"-i-|-l 

Then multiplying (|20|l by z'^ and summing over d > 2" 
yields 



d 



[G{z,r)-H{z,T)] 

oo 2" 

= E E Pid-d',r) 

d=2" + l d'=2"-i + l 
2" oo 

= E E Pid",r) 

d"=2"-i + l d'=2'>-i+l 

= -G{z,t)^H{z,t) 



d_ 
97 



P(d',r) 



(23) 



2"-i<ci'<2" 



where we have used that in stage n of the dynamics 
P(d",T) = for d" < 2"-i. We may integrate ^ 
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FIG. 2; Domain length distributions Pn{d) at the end of stage 
n — 1 of the low T coarsening dynamics, for initial temper- 
ature T — oo. Open symbols and lines; Theoretical results, 
calculated from 1271 1. for n — (Oi initial condition), 1 (□), 
2 (0), 3 (A). Full symbols; Simulation results for a chain 
of length L = 2^^ and e = 10"* {n = 1, 2) and e = 10"^ 
(n = 3). Inset; Scaled predictions 2"~^P„(d = 2"~^x) vs. x 
for n = 1, . . . , 8. Bold line; Predicted scaling function 13H . 



by rewriting it as 
1 



d 



[1-G(z,t)] dr 



G{z,t) 



d 

OT 



giving 



1 - G(z,0) 



exp r) + 0)] 



(24) 



(25) 



At the end of stage n, corresponding to r — > oo, all 
domains that were active during that stage have disap- 
peared, and so H{z, oo) = 0. Thus 



G{z, oo)-l = [G{z, 0) - 1] cxp[H{z, 0)] 



(26) 



Now recall that we are considering stage n of the dynam- 
ics. The initial condition for stage n + 1 of the dynam- 
ics will be given by the distribution P(d,t) at the end 
of stage n. Thus defining G'„(z) = G(z,0) for stage n, 
with a similar definition for the active generating func- 
tion Hn{z), Eq. H2t)|) relates the different stages of the 
dynamics through 



G„+i(^) - 1 = [G„(z) ~ 1] exp[i7„(z)] 



(27) 



This exact result connects, through their generating 
functions, P„(d) and Pn+i{d) which are defined as the 
domain length distributions at the end of stages n — 1 
and n of the dynamics respectively. It can be checked 
that for n — 0, where the system does not cross energy 
barriers, Eq. H27() gives results compatible with the known 
exact solution of the T = dynamics . 

Iterating (|27|) from a given initial distribution Paid) 
gives Pn{d) for all n = 1,2,... Fig. |21 shows numerical 



results for the case where Po{d) is the equilibrium dis- 
tribution (jSJ corresponding to an initial temperature of 
T = oo. It is clear that a scaling limit emerges for large 
n, i.e. that the rescaled distributions 



P„(2;) = 2""ip„(d) where x 



(28) 



converge to a limiting distribution P{x) for the scaled 
domain size x. This is just a statement of the invariance 
of the coarsening processes in each stage once the domain 
sizes are rescaled by the characteristic size 2"~^. 

The change to a continuous variable x for the do- 
main lengths simply results in the generating functions 
Gn{z), Hn{z) being replaced by Laplace transforms 5n(s) 
and hn{s), via the correspondence z'^ e~*^. Invariance 
of the scaled domain distribution under H27|) then gives 
the equation 



g{2s) - 1 - [gis) - 1] exp[/i(s)] 



(29) 



where 
9is) 



dxP{x)e-'^, his) 



dxP{x) e 



We found a solution to H29(l by noting that the numerics 
strongly suggest P{x) = 1/a: for 1 < a; < 2. Using this 
as an ansatz implies 

h{s) = Ei(s) - Ei(2s) where Ei(s) = / du 

Js U 

which when inserted into 1)29(1 yields 

[1 — g(s)] exp[Ei(s)] = const. 

The requirement that 17(5) ^0 for large s fixes the con- 
stant as unity, giving 



g{s) = 1 - exp[-Ei(s)] 



(30) 



Expanding the exponential as a series allows the Laplace 
transform to be inverted term by term and one obtains 



74-200 



dse^^g{s) 



7 — zoo 
-j^ ^7+joo 

27ri 



\m+l 



dse'' 



m 



- (_l)™+i .00 ™ ^ / ™ N 



\k=l 



e(a; - 1) - - e(a; - 2) — ^ '- 

X x 



(31) 



where Q{x) is the Heaviside step function. This result 
for the scaling function is shown in Fig. |21 and shows 
good agreement with the numerically calculated Pn{d) 
for large n. Note that the series (|3If) has singularities in 
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the fc-th derivative at the integer values x = k + 1, k + 2, 
. . . , a fact whose physical origin we discuss further in 
the next subsection. The average domain length in the 
scaling limit is given by dn = 2'^~^x; from the results for 
P{x) we find x = exp(7) = 1.78 . . ., where 7 is Euler's 
constant. 



D. A family of related coarsening models 

Surprisingly, the scaling function H31() for coarsening 
dynamics in the East model is identical to that for the 
'paste-all' model of coarsening dynamics, where at each 
step in the dynamics the shortest domain is selected and 
'pasted' onto its left or right neighbour [s^l- To under- 
stand this, we decide for definiteness that domains are to 
be pasted to the right. Then the paste-all model can be 
obtained from a modification of the hierarchical coarsen- 
ing dynamics discussed so far — one merely needs to as- 
sume that the rates T{d) are now all well separated from 
each other, so that at any given stage in the dynamics 
only domains of a single size d are active. It is then natu- 
ral to consider a family of models which interpolates be- 
tween paste-all and the East model, by assuming that the 
active domain sizes d at any stage n — whose coalescence 
rates r(d) are comparable to each other but well sepa- 
rated from all other rates — are those with a"~^ < d < a", 
with a a constant in the range 1 < a < 2. 

All of the arguments in the previous section apply to 
this modified model if powers of 2 are replaced by powers 
of a in the appropriate places. In particular, Eq. (|29|) for 
the Laplace transform of the scaling distribution becomes 

5(as)-l = [5(s)-l]exp[;ia(s)] (32) 

with 

ha{s)= dxP{x)e-''' (33) 

generalizing the earlier h{s) = h2{s). We now show 
that H32|) has the same family of solutions for any a < 2, 
thus proving in particular the link between the scaling 
distributions of the East model (a — 2) and the paste-all 
model (a 1). To see this, consider first a ^ 1; we 
should then recover the scaling equation for the paste-all 
model. Indeed, setting a = I + S and expanding to first 
order in d, Eq. (|32|l becomes 

sg'is) = [g{s) - l]P(l)e-^ (34) 

which is equivalent to Eq. (38) in Ref. 32] . Integrat- 
ing gives 32] gis) = 1 - exp[-P(l)Ei(s)], and the 
Laplace transform can be inverted as in (|31|l . giving 

p{x) = 6(2: - 1) — - e(x - 2) ^-^ i '- + ... 

X X 

Given that only the first term contributes for a; < 2, it is 
now easy to verify that this is also a solution of H32|) for 



any a < 2: from H33|l one has ha{s) = — P(l)[Ei(a,s) — 
Ei(s)] =ln[l-g(as)]-ln[l-5(s)]. 

In summary, we have shown that there is a whole fam- 
ily of coarsening models which interpolate between the 
paste-all and East models and which have the same scal- 
ing solution. The solution is, in principle, parameterized 
by P{1) but, as discussed in ^z], the requirement of a 
finite mean for the scaling distribution imposes P(l) — 1. 
Thus g{s) = 1 — exp[— Ei(s)] as we found in (|30|l . 

Physically, the common feature of all the models in the 
family is that active domains can never be created during 
the stage where they are active; this breaks down for 
a > 2 where indeed the scaling solution no longer applies. 
This argument clarifies the origin of the singularities at 
integer arguments of P{x): the existence of a shortest 
scaled domain length x = 1 implies that the shortest 
inactive domain that can be created has x — 2^ and this 
effect then propagates to x = 3, 4, . . . as the dynamics is 
iterated. We initially thought that the fact that in the 
East model the active domains cover a range x £ [1,2] 
bounded by two integers played a role, but this is clearly 
not so since for general a < 2 this range becomes x G 
[l,a]. 



E. Anomalous coarsening 

Having solved the dynamics in terms of the different 
stages labelled by n, we now translate these results into 
actual time dependencies. As before, we consider a sys- 
tem quenched from, say, T = oo, to a temperature T 
at time i = 0. If data for e.g. the average domain length 
d are plotted against the scaled time variable ly — Tint, 
then for T the n-th stage of the dynamics shrinks 
to the point v = n. In this limit we predict that, for 
n — 1 < v < n, the domain length distribution is Pn{d) 
as defined by the recursion (|27|) . The average domain 
length div) will follow a 'staircase' function, jumping at 
V = n from c?„ = ^j^Pnid)d to dn+i. This was illus- 
trated and verified by low-T simulations in |18|. 

For large n, i.e. in the scaling regime of large d (but still 
d <C dcq) we know that dn = 2"-~^x where x = 1.78 .. . 
The staircase function d(i/) is therefore bounded between 
2''-'^x < d < 2''x, giving i < J/(xF'"2) < i ^hen ex- 
pressed in terms of ordinary time t. This shows that 
the hierarchical dynamics gives rise to anomalous coars- 
ening, i.e. the typical domain size grows as d ~ f^'"^, 
more slowly in time than a usual power law with T- 
independent exponent. The exponent can also be de- 
duced from the scaling t ~ r(c?)^^ ^ e" = exp(rt/T) 
of the time required to eliminate a domain of the typical 
size d. Since d ~ 2" asymptotically, one has n « Ind/ ln2 
and this time scales as t ~ d}^^^'^'^, giving J ~ t"^'"^^ as 
before. For decreasing T the coarsening becomes anoma- 
lously slow and in fact logarithmic for T — > 0. 

By extrapolating the anomalous coarsening law to the 
equilibrium domain length d^q = exp(l/T)-t-C'(l), we can 
estimate the equilibration time of the system for T — > 
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as teq - = exp(l/T2ln2). Of course, in performing 
this extrapolation, we are using the anomalous coarsen- 
ing law in the regime where domain sizes are 0{l/e) and 
the assumption of irreversible down- flips no longer holds. 
Given this we could have equally well extrapolated d ^ 
t"^'"^ to some multiple adcq of the equilibrium domain 
length, which would yield foq = exp(l/r2ln2 - A/T) 
with A = — In a/ In 2. With the sign as defined, A should 
be positive, corresponding to a < 1; a negative value of A 
is excluded as it would violate the relaxation time bound 
of jl^ . Thus we expect to give only the dominant di- 
vergence of the equilibration time teq, with subdominant 
correction factors exp(— A/T) making teq <C i* generi- 
cally. Our analysis in Sec. lIVI suergests. on the other hand, 
that the longest relaxation timescales at equilibrium are 
given directly by . This is not unreasonable: the time to 
reach equilibrium from an initial high temperature con- 
figuration is expected to be much shorter than the time 
for correlations to be erased starting from an equilibrium 
(low temperature) initial condition, simply because the 
typical domains to be eliminated are shorter. 

Finally, we comment briefly on the implications of the 
coarsening dynamics for the two-time spin autocorrela- 
tion function C(t, iw) = (si(t)sj(tw)) - (sj(t)) (sj(iw)) 
where t > t„. Since in the limit T — > spins flip down 
irreversibly, Si(t) = 1 implies that Si(iw) = 1 at the ear- 
lier time tw, thus C{t,t^) = c{t)[l — c{t„)] where c{t) 
is the time-dependent up-spin concentration. As a func- 
tion of t, C(t,t^) will thus exhibit the same plateaux as 
c(t\ and these were indeed observed in the simulations 
of [23 • Also, the normalised autocorrelation becomes 
simply C(t,tw)/C(iw,tw) = c{t)/c{t^) = d{t^)/d{t). 
This dependence on only the ratio of the relevant length- 
scales at times and t is natural for a coarsening pro- 
cess, and gives a good fit to the data of p^ . 

IV. EQUILIBRIUM DYNAMICS 

A. Functions of interest 

We now turn to the equilibrium dynamics of the East 
model at low temperatures. These are encoded in re- 
laxation functions such as correlations and spin persis- 
tence (see below); response functions provide no new in- 
formation since they are related to correlations via the 
fluctuation-dissipation theorem (FDT). For correlation 
functions standard choices would be the local spin au- 
tocorrelation {si{t)si{0)) — (? or the correlation function 
of the up-spin concentration (or magnetization) M(t) = 
(1/L) Si(t). However, the directionality of the kinetic 
constraint can be shown to imply that all non-local cor- 
relations (s,;(t)sj(0)) ~c? vanish |33| so these two choices 
give the same information. We can thus focus on the 
normalized autocorrelation function 

^ ^ (g,(Og,(0))-C^ ^ P(g,(t) = l|g,;(0) = l)-C 

c — (? 1 — c 

(35) 



This decays from C(0) = 1 to C(t ^ 00) = 0; the second 
equality in l|35|l implies that it is essentially the condi- 
tional probability for an up-spin that was up at time 
also to be up at time t. Closely related is the up-spin per- 
sistence function T'i{t), which gives the probability that 
a spin is up throughout the time interval [0,t], and the 
analogously deflned down-spin persistence Voit). 

Fortunately, it turns out that C, Vi and Vq do not 
need to be analysed separately because they all become 
identical in the limit T — > (taken at a constant value 
of, say, C{t)) that we are interested in. To see this, note 
that since in the East model the spins are only coupled 
via the kinetic constraint, which acts to the right, the 
dynamics of a spin Si cannot influence that of its left 
neighbour Si_i. (This is strictly true only in the limit of 
an inflnite chain L — > 00, but that is precisely the case of 
interest.) Furthermore, the dynamical evolution of spin 
Si is that of a single spin in a field whenever Si_i — 1, 
and is completely frozen otherwise. It follows that if, 
between times and t, Si-i has been up a fraction m of 
the time, then the evolution of Si has been that of a single 
spin over time mt. But the single-spin correlation and 
persistence functions are trivial to work out; denoting 
the distribution of m for a given t by P{m]t), one thus 
has the simple expressions 

C{t) = f dmP(m;t)e-^^+''>"'* (36) 
Jo 

Vi[t) = [ dmP{m;t)e-'^* (37) 
Jo 

Vo{t) = [ dmP{m;t)e-'"'* (38) 

It follows trivially that 

Cit)<Pi{t)<ro{t) (39) 

Now since P{s,(t) ^ l\s,{0) ^ 1) > Vi{t) (the probabil- 
ity for a spin to be up at times and t must be greater 
than that for it to be up at and t and all intermediate 
times), one has from (I35|l 

{l-c)C{t)+c>Vi(t) (40) 

Eqs. (|39I4U|) together show that C(t) and Vi{t) coincide 
in the low-T limit (c w e — > 0) whenever the values of the 
fmictions themselves are of 0{1). 

To get the analogous result for C and Pq, consider 
P{m;t). The concentration of persistent up-spins and 
down-spins in the system is cPi(t) and (1 — c)PQ{t), re- 
spectively, and these have m = 1 and m — 0. Thus 

P{m;t) = {l-c)Po{t)d{m) + cPi{t)S{m-l) + ... (41) 

where the dots indicate contributions for m strictly be- 
tween and 1. Inserting into (|37|) and separating off the 
(5(r7i)-term gives 

p,{t) = {l-c)Po{t)+ f dmP{m;t)e^"'' > {l-c)Po{t) 
J>o 
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FIG. 3: Simulated MFTPs r^fp for down-flips (solid) and r^fp 
for up-flips (dashed) as a function of unsealed distance d, for 
e = 0.3,0.2,0.1,0.05,0.02 (bottom to top). The thin solid 
lines have slope 1, showing that for large d the MFPTs vary 
linearly with d. Error bars are too small to show on the scale 
of the plot (< 3.5% relative error on Tmfp and r^fp). 



Together with ^ it follows that also Pi{t) and Po{t) 
coincide for c — > 0. Thus C = Vi = Vq in the limit, and 
we can restrict attention to e.g. the up-spin persistence 
Vi{t) in the following. 

We stress once more that the directionality of the ki- 
netic constraint is essential for the above arguments to 
work. In the undirected version of the East model, i.e. 
the one-dimensional Fredrickson-Andersen model with 
/ = 1, it would still be true that the dynamics of spin 
Si is determined by the amount of time which Si-i (and 
Si+i) spend in the up-state. But this time cannot be de- 
termined independently of the state of spin Si, since Si 
itself affects the dynamics of Si-i (and st+i). 



B. First passage times 

To get some insight into the equilibrium dynamics, we 
generalise the domain coalescence rates r(d), discussed 
in Sec. IIIII for d = 0(1) and e — *■ 0, to the regime of 
domain sizes typically found in equilibrium, d ~ dcq ~ 
1/e. Starting from the state (sq, Sij • ■ • j sj) = 10 . . . 01, 
we define a mean first passage time (MFPT) Tjnip{d) as 
the mean time for the spin Sd to fiip down for the first 
time, and we set T{d) = T"~fp(d). Note that in Sec. lIIIBl 
we required for a 'first passage' that not just Sd but all 
other spins si, . . . , Sd-i be down; here we just require 
that Sd — 0. In the regime considered in Sec. IIIII the 
two definitions are equivalent to leading order, since the 
relaxation process has passed its 'highest point' once Sd 
has flipped down, resulting in the same energy barrier 
whether or not the process continues to the empty state. 

Consider now the dependence of rnifp(rf) on d. From 



the discussion of the out-of-equilibrium dynamics we 
know that T^{p{d) exhibits a step-like variation with d 
for d = 0(1) and small e, varying only by factors of order 
unity within each range 2"~^ < d < 2" but increasing by 
~ between ranges. As n increases at fixed e, however, 
the 0(1) changes within ranges become larger and even- 
tually comparable to e~^; T^nipid) will then increase more 
smoothly with d. The numerically simulated MFPTs in 
Fig. 121 confirm this. Finally, for d larger than deq ~ 1/e, 
the figure shows that the MFPTs increase linearly with 
d. Intuitively, one then has a 'front' of up-spins that 
creeps 'ballistically' along the chain. This can be moti- 
vated by arguing that, for d above the typical equilibrium 
lengthscale deq, the spin front no longer remembers that 
it originated from a single up-spin a distance d to the 
left. Instead it behaves as on an infinite chain, where its 
propagation velocity is necessarily constant. 

So far we have considered the variation of rmtp (d) with 
d at fixed e. One suspects, however, that e.g. the decay 
of the persistence function Vi {t) is governed by a length- 
scale d{t), giving the size of the largest domains that have 
had time to equilibrate, and more precisely — since we are 
considering equilibrium dynamics — by the ratio d{t)/deq. 
We thus define scaled domain sizes d — ed {= d/dcq for 
e — > 0) and now ask how rmtp (c?) scales with T for fixed d. 
For small d, extrapolating from the regime of d — 0(1), 
and using n « In d/ In 2, one expects 



-lnd/ln2 _ ^- In d/ In 2- 1/(T In 2) 

/ Ind 1 



exp 



rin2 T^Iti2 



Rearranging, one has 



Tmfp(rf) 



Tin 2 



From the downward curvature of the lines in Fig. 13 one 
sees that this expression will not be valid for larger d, 
where the r.h.s. will cross over to a slower increase with 
d. We thus make a time scaling hypothesis which replaces 
d by a more general scaling function /(d) of d, i.e. we 
assume that the rescaled MFPT approaches 



- Tmfp(d) 



Tin 2 



/(d) for r ^ (42) 



when the limit is taken at fixed d. Fig.QJleft) represents 
numerical data for the MFPTs in this scaled form and 
shows that the assumption 1)42(1 is certainly plausible for 
not-too-large d. For larger d, convergence to a scaling 
limit is not yet evident from the data that we can gen- 
erate on practical simulation timescales. This is because 
for d ^ doq the ballistic propagation discussed above im- 
plies fijifp oc d"'" ''^ ^ (X d^ ^ which will tend to a constant 
as T — > but does so very slowly. For the scaling func- 
tion /(d) this implies that it must approach a constant 
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FIG. 4: Left: Logarithmic plot of scaled FPTs f, as defined 
in II42II . versus scaled distance d, for e = 0.2,0.1,0.05,0.02 
(top to bottom on right). The dotted line shows a qualitative 
sketch of the scaling function f{d), with / — » 1 (In / — » 0) for 
d — > oo; see text. Right: Up-flip (dashed) versus down-fiip 
(solid) FPTs and their fluctuations, versus e for constant d = 
ed — 0.1. Thin lines and error bars: mean and ±i(standard 
deviation) of the log-rescaled FPT Inf. Thick lines: Infmfp. 
As e decreases, the fiuctuations in Inf are seen to decrease: 
error bars shrink, and Infmfp and the mean of Inf become 
closer. Also, up-flip and down-flip times become closer; both 
observations support our scaling hypothesis for the FPTs. 



correct intuition is that both up- and down-spins only 
become mobile once an up-spin 'front' from an up-spin 
to the left has reached them, and that for T ^ the time 
required for this front propagation vastly dominates the 
effect of the different flip rates of up- and down-spins once 
the front has reached them. Fig. 01^ right) again supports 
this assumption with numerical data for d = 0.1. 

In summary, we assume in the following that (|42|l is 
valid for unaveraged up-flip and down-flip FPTs, with 
the same scaling function f{d). A consequence of this 
is continuous timescale separation: for any two rescaled 
distances c?i < d2, the ratio of the corresponding FPTs is 
[/(d2)//(di)]^/^^'"^) and diverges as T ^ 0. In the limit, 
this means that the equilibration of domains of length di 
proceeds infinitely more quickly than for any even slightly 
larger length c?2. This insight is key for the construction 
of the superdomain model described next. A proviso is 
that we have assumed here that the scaling function f{d) 
is strictly monotonically increasing with d. In principle, 
it is possible that f{d) could increase monotonically only 
up to some finite d*, and be exactly constant thereafter. 
This seems to us implausible, however — e.g. it is difficult 
to conceive of a physical mechanism causing the singular- 
ity at — and, as discussed below, would also give very 
unusual predictions for the time-dependence of the low-T 
relaxation functions. 



as d — > OO] recall that this limit is taken after the limit 
T ^ at fixed d. We note that the bounds of 12] imply 
that the MFPTs rmfp cannot exceed t^, by more than fac- 
tors of 0(1); this imphes from |jl21l that f{d) < 1. Our 
numerical data in Fig. ^Jleft) are restricted to values of 
T which are still too large to determine f{d) with any 
accuracy, but the tentative extrapolation to T — > indi- 
cated in the figure is certainly consistent with the bound 
f{d) < 1. 

We work below with a slightly strengthened version 
of (|42|l . The first passage-time t for down-flips is a 
fluctuating quantity, with mean r^fp. We assume that 
fluctuations in t are small enough so that 142|) holds 
even for the fluctuating r, which means that for T — > 
the rescaled FPT f{d) = (r(d)/i*)'^'"2 becomes non- 
fluctuating. This assumption is not as strong as it may 
sound; because of the exponentiation by Tin 2 it holds 
e.g. if the relative fluctuations of the unsealed FPT 
t(c?) remain of 0(1) as T ^ 0. Further confirmation 
comes from Fig.QJright) which shows numerical data for 
d = 0.1. 

In a final generalization we also assume that H42|) ap- 
plies if we consider FPTs t' for up-flips, where we start 
with the empty state 10 ... 00 and r' is the first time 
where the last spin, Sd, has flipped up. Naively one 
might have suspected that t' and r differ by a factor 
~ 1/e. However, the equivalence oi'Pi{t) and Po{t) that 
we proved for low T already suggests otherwise. The 



C. Superdomain model 

We now exploit the idea of continuous timescale sepa- 
ration to construct an effective description for the low-T 
equilibrium dynamics of the East model. It is natural 
to work with the equilibrated lengthscale that corre- 
sponds to the relaxation timescale t,- we are considering; 
from (|42|l . the two are related by 



Here and in the following, the limit T —> (or equiv- 
alently e — > 0, or c — > 0) is always understood. For 
simplicity we will simply call c?r 'time' where there is no 
ambiguity. We stress that we work with rescaled length- 
scales d — de throughout; a value of d of over unity thus 
corresponds to a very large domain length d = d/e in the 
limit e — > 0. 

Consider now an equilibrium configuration of the spin 
chain, consisting of up-spins separated by long domains 
of down-spins. From the (rescaled) domain sizes are 
distributed according to a simple exponential, P{d) — 
exp(— d). A 'time' d^ later, each of the up-spins will 
have developed an equilibrated zone of length dr to its 
right: from continuous timescale separation, the FPT to 
reach any spin within this zone is much shorter than t^-. 
Throughout the zone we should thus indeed have equili- 
brated spins, which are independently up with probabil- 
ity c = e/(l + e) ~ £. The domains defined by up-spins 
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within these equihbrated zones will not be of interest in 
the following. Instead we focus on the 'superdomains' 
which are naturally defined by the unequilibrated up- 
spins, to which we refer as 'superspins'. Thus each su- 
perdomain is bounded on the left by such a superspin, 
followed by an equilibrated zone of length di- and then a 
string of down-spins. 

To understand the dynamics of superdomains as we 
look at increasing times d^, it is easiest to consider a sys- 
tem for which the rescaled chain length eL and therefore 
the typical number of superdomains is large but finite; 
there is then at all times a nonzero minimum superdo- 
main length. On increasing d,- from zero, all superspins 
remain as they are as long as the equilibrated zone of 
each superdomain has not yet reached the superspin to 
the right. However, when d^ becomes equal to the small- 
est superdomain length in the system, di say, then this 
superdomain's equilibrated zone 'catches up' with super- 
spin 2 bounding the next superdomain (of length ^2) on 
the right. This superspin now becomes equilibrated, and 
so one might assume that the two superdomains just co- 
alesce, forming a single superdomain of length di + d2. 
This, however, will only be the case if there are no up- 
spins in the equilibrated zone of superspin 2 at this 
time; since this zone has length dr, the probability for 
this event is exp(— dj.). Otherwise, i.e. with probability 
1— exp(— dr), there will be at least one up-spin in the equi- 
librated zone of superspin 2. The leftmost of these, hav- 
ing lost the superspin from which it became equilibrated, 
becomes frozen and thus itself turns into a superspin, 2'. 
The distance S between 2' and the old superspin 2 has 
probability distribution exp(— (5)/[l — exp(— dp)] over the 
interval S G [0, dj], and the two new superdomains have 
length d[ — di + 6 and d'2 — d2 ~ 5 (see Fig. [S)). 

What makes the superdomain model nontrivial is that 
in the second case the process of superspin elimination 
and regeneration can now continue. Superspin 2' will 
immediately (on the timescale dr being considered) equi- 
librate a zone of length d^-. If d^ < d^, then this zone 
includes superspin 3, which now becomes equilibrated 
along with a segment of length d^ — d^ to its right. This, 
however, still leaves a segment of length d^ — {d^—d'2) = d^ 
of the old equilibrated zone of superspin 3, which is now 
frozen (see Fig. EJ. If there is no up-spin within this 
segment when superspin 3 equilibrates, then superdo- 
mains 2' and 3 coalesce; this happens with probability 
exp(— ^2). Otherwise, the leftmost up-spin in the seg- 
ment freezes into a new superspin 3', and the process 
continues by iteration. It is clear that only superdomains 
of length d > dr have been generated when the process 
terminates. One can thus now increase dr until the new 
minimum superdomain size is reached, at which point a 
new relaxation process as described above starts. 

Above, we arrived at the superdomain model start- 
ing from the hypothesis of continuous timescale separa- 
tion. The model lets us access, via an effective descrip- 
tion, very low temperatures corresponding to extremely 
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FIG. 5: Illustration of superdomain dynamics. Superspins 
are shown by the thin vertical rectangles; each has an equili- 
brated zone of length dr to its right, indicated by a horizontal 
rectangle. Shown is a situation where the relaxation 'time' 
dr has just become equal to the smallest superdomain length, 
di. The equilibrated zone of superspin 1 then 'catches up' 
with 2 and equilibrates this spin (first line). At this moment, 
there can be up-spins within the former equilibrated zone of 
2; if there are, the leftmost of these becomes a new super- 
spin 2', bounding a superdomain of length d'2 (second fine). 
If d'2 < dr, the newly created equilibrated zone of 2' catcfies 
superspin 3. Within the zone of length d'2 that is now no 
longer equilibrated (dashed) up-spins can again remain, with 
the leftmost becoming a new superspin 3'. In the example, 
the relaxation process at dr stops at this point since d'^ > d^; 
otherwise, it would continue in the same fashion. 

long relaxation times. By definition it is therefore difh- 
cult to demonstrate superdomain-type dynamics on the 
timescales of a numerical simulation. Nevertheless, the 
sample run at e = 0.02 shown in Fig. illustrates some 
important features of superdomain dynamics, in particu- 
lar the regeneration of superspins and the resulting prop- 
agation of equilibrated zones. 

To summarize the superdomain model, let us restate 
its dynamics. We present this in the form of a schematic 
simulation algorithm; a formal definition of the stochas- 
tic evolution of the sequence of superdomain lengths is 
of course possible but would be more awkward. We re- 
emphasize that all lengths are rescaled lengths, d = de. 

1 . Initialize superdomain lengths from an exponential 
distribution P{d) = cxp(— d). 

2. Set dr — size of smallest current superdomain. Let 
that size be di, with d2, ds, . . . the sizes of the su- 
perdomains on the right. Set i — 1. 

3. Delete the now equilibrated superspin i + 1. With 
probability exp(— di), the superdomains i and i + 
1 coalesce and the relaxation process at this dr is 
complete: let di ^ d^ -f di+i, delete superdomain 
i + 1 and go back to step 2. 

4. Otherwise, a new superspin (i + 1)' is created, a 
distance S to the right of the old one which is dis- 
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FIG. 6: A direct simulation run on a system at e = 0.02 which provides support for the notion of superdomain-type dynamics. 
Shown is a section of 180 spins out of a much longer chain, for a single simulation run. The lines correspond to seven evenly 
(logarithmically) spaced time intervals within the range 6.26 X 10** < t < 10^". The boxes indicate the current magnetizations 
Si of the spins, determined by averaging over the relevant time interval. A logarithmic scale is used, so that the highest boxes 
correspond to Si = 1, the dashed lines to the equilibrium value Si = c, and the baselines to Si — (? . (Lower values of Si are 
not shown.) Filled boxes indicate spins which are persistent over the time interval (si = 1). Note the similarity between the 
first two lines here and the superdomain dynamics sketched in Fig. |5] within the time interval (ii), the equilibrated zone of 
superspin 1 catches up with 2; a new superspin (2') is created and eliminates superspin 3. Lines (v) and (vii) show events where 
superspins 1 and 4 are eliminated. Of course, e here is still too large to be in the asymptotic limit T — ^ where superdomain 
dynamics applies exactly. This is why, in contrast to Fig. |S] the equilibrated zones of the superspins are not all of the same 
length and do not have sharply defined boundaries. 



tributed according to P{5) oc exp(— (5) over the 
range (S G — di,dr]. This gives new superdo- 
main lengths: <— + (5 and d^+i <— d^+i — (5. If 
dj+i < dr, increase i <— z + 1 and go back to step 
3. Otherwise the relaxation process at this is 
complete; go back to step 2. 

Ideally, one would like to solve the above superdomain 
model directly and thus compute the dependence of the 
persistence function on dr analytically. The iterative pro- 
cess of superspin elimination and regeneration is difficult 
to keep track of, however, and we have not been able to 
find an analytical solution. Nevertheless, because the su- 
perdomain model expresses times in terms of the length- 
scales dr, which unlike the timescales themselves do not 
diverge for T — > 0, it can be simulated easily and ac- 
curately; we have used systems of 10^ superdomains (at 
dr = 0) , checking that finite-size effects are negligible and 
typically averaging our results over 10^ simulation runs. 

The key quantity that we want to predict from the 
superdomain model is the persistence function. The up- 
spin persistence V\ is the fraction of up-spins that have 
never flipped since dr — 0. Since at dr = all up-spins are 
superspins, V\ is the fraction of these initial superspins 
that have never been equilibrated. We show the results in 
Fig.d V\ initially decreases linearly with dr, but then the 
decay becomes much faster and indeed super-exponential 
in dr. 

If, as we have claimed, the superdomain model is the 
correct effective description for the T — > dynamics of 
the East model, then it must obey the exact identity 
V\ = ■Po derived in Sec. IIV Al Vq is measured in the su- 
perdomain model as the fraction of the chain which has 




FIG. 7: Predictions of the superdomain model. Main plot, 
dashed lines: Response functions for eo/e = 0.95 and eo/e = 
1.05 versus relaxation 'time' dr; on the latter, a few error bars 
are shown where they are significant. Solid lines: Persistence 
functions Vi and Vo\ the two lines are indistinguishable by 
eye. The inset shows the scaling function g{t) for the equilib- 
rium relaxation functions, as derived from the superdomain 
persistence and the scaling function for the FPTs sketched in 
Fig.mieft). 



never been swept by an equilibrated zone, a quantity that 
one might not have naively suspected to be connected to 
the number of persistent superspins. Nevertheless, our 
simulations show that indeed Vi = Vq to very high ac- 
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curacy (~ 1%, less than the relative error on the mea- 
surements of Vi and Vq) in the superdomain model. In 
fact, both quantities are plotted in Fig. [3 but are indis- 
tinguishable by eye. This provides strong support for the 
correctness of the superdomain model; if e.g. the mech- 
anism for regenerating superdomains is neglected, one 
finds that the condition Vq = Vi is violated. 

A further consistency check on the model is obtained 
from the requirement that it should represent the equi- 
librium dynamics. The concentration of up-spins in the 
system, i.e. the magnetization M = i^/L)'^^Si, should 
thus remain independent of dy. Within the superdomain 
model, if there are Ns superspins at time dr, one has 

M = N,il + cd,/e)/L « iNjL)il + d,) (43) 

since in addition to the superspins there are a further 
c X Ngdr/e up-spins on the Ngdr/e sites covered by the 
equilibrated zones. This quantity should equal M = c « 
e independently of d^. Hence the average (rescaled) dis- 
tance between superspins, eL/Ng, which is identical to 
the average superdomain length, must equal 1 + d^ at 
time dy. This is indeed what our simulations of the su- 
perdomain model show. One can push this comparison 
further and consider not just the total up-spin concentra- 
tion, but the distribution of sizes of the domains formed 
by all spins (rather than just the superspins). This distri- 
bution must be a pure exponential, independently of di-. 
It can be expressed in terms of the superdomain distri- 
bution by an appropriate convolution which accounts for 
the fact that additional up-spins exist in the equilibrated 
zones. Omitting the details, we only state that one finds 
in this way that the superdomain size distribution at time 
dr must be 

P(d; d~.) = e(d - d,)c-^^-^'^ (44) 

Again we find that this is verified in our simulations. 
Note that l|44|) has a simple intuitive interpretation: it 
corresponds to an exponential distribution of the seg- 
ments of down-spins separating the equilibrated zone of 
each superdomain from the next superspin on the right. 

A final check on the superdomain model is that it 
should obey FDT — since the original model obeys de- 
tailed balance FDT is automatically satisfied but this is 
not guaranteed for the superdomain description. As ex- 
plained in Sec. lIV A1 all non-local correlation (and hence 
response) functions vanish in the East model, so that 
one is free to consider either a local response of Si to 
a local field, or a response of the magnetization M to 
a uniform field h. Choosing the second option, the en- 
ergy function is modified to E' = (1 — h) Si which 
is equivalent to changing e to e' = eexp(/i/T) (or tem- 
perature from T to T / {1 — h)). The response to a field 
h switched on in the distant past and switched off at 
t — Q can thus be measured by initializing the system 
in an equilibrium state corresponding to e' and mon- 
itoring the evolution of M during the subsequent dy- 
namics at e. By FDT this switch-off response should 



have the same time-dependence as the correlation (and 
hence the persistence) functions. In the superdomain 
model, the measurement is performed by initializing the 
superdomains with a modified domain size distribution 
P(d) = (e'/e) exp[— (e'/e)d] and then tracking the de- 
cay of the magnetization, measured as in (|43|l . to its 
equilibrium value. The response functions simulated for 
e'/e = 0.95 and 1.05 are plotted in Fig. [7| above, and 
show that the superdomain model indeed obeys FDT. 

Pushing the above scenario further, one could consider 
nonlinear responses in the superdomain model, in partic- 
ular a large ratio e'/e which corresponds to a quench to a 
much lower temperature. (For the superdomain model to 
remain applicable we still need e' ^ 1, of course.) The 
initial scaled domain lengths in the system are then of 
order e/e' ^ 1. In the regime of small dr where these do- 
main lengths are removed from the system, one sees that 
the probability 1 — exp(— dj-) ~ dr of creating a new su- 
perspin is very small, tending to zero for dj. e/e' 0. 
In this limit the superdomain dynamics becomes 'take 
the smallest superdomain and coalesce it with its neigh- 
bour on the right', which is precisely the paste-all model 
discussed in Sec. IIII DI While superdomain sizes remain 
^ 1, the domain size distribution will thus be driven to 
the scaling distribution of the paste-all model. As demon- 
strated in Sec. IIII D1 this scaling distribution is the same 
as that for the coarsening dynamics of the East model. 
Thus, the prediction of the superdomain model for the 
form of the domain size distribution, a long time after a 
quench from e' to e with 1 3> e' 3> e, matches up pre- 
cisely with that predicted in Sec. IIII CI for a quench from 
1 ~ e' e. (Note that in the first case we are in princi- 
ple talking about the superdomain size distribution, not 
the actual domain distribution as in Sec. IIII (Jl However, 
for small di- the two are identical since the number of up- 
spins within the equilibrated zones of the superdomains 
are negligible.) 

Writing our result for the persistence function in the 
superdomain model (either Vi or Po, since they are iden- 
tical) as Psd(dr), we can now translate it into a scaling 
prediction for the time-dependence of the correlation and 
persistence functions at low T in the East model. Using 
the inverse of the scaling function /(d) from (|42|l . one 
has 

cit) = v,{t) = r,{t) = irHi)) , f = (^-L^ 

(45) 

Thus the relaxation functions show strong stretching, 
with the stretching exponent T In 2 decreasing to zero for 
T ^ 0. However, the scaling function g{i) — 'Psd(/~^(i)) 
is non-exponential. In fact, since /(d) approaches a con- 
stant /oo for d — > oo as argued after the scaling 
function in H45|l decays to zero at a finite value t = foe ■ 
On the basis of the results ofp^ we would conjecture this 
value to be /oo = 1, see after (|42|l : the scaling function 
g{t) for this case is sketched in the inset of Fig. [7| It 
is important to bear in mind that the limit of T ^ 
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FIG. 8: Comparison of up-spin persistence function Vi from 
simulations with prediction of superdomain model, for e = 
0.02. A log- log plot of — InPi(t) is shown, which would be 
straight for stretched exponential relaxation. The inset shows 
directly. Thick lines: Vi as simulated in three runs for 
a chain of length L = 2^"^; the correlation function C{t) is also 
shown for one run but is indistinguishable from Vi {t) as pre- 
dicted for low e. Thin lines: Predictions of the superdomain 
model (see text). Dashed/solid: t determined from FPTs for 
up- and down-flips, respectively; left/right curve of each pair: 
using mean of the log-FPT and the log of the MFPT, respec- 
tively. 



considered here is taken at constant or, equivalently, 
constant t. If we look instead at fixed nonzero T then the 
relaxation functions are of course nonzero for all t, but 
this does not contradict g{t) dropping to zero at t = foo- 
To see this, note that the asymptotic decay of all relax- 
ation functions should be to leading order an exponen- 
tial with the longest relaxation time, ^ exp(— This 
does remain nonzero for all t; but in terms of the scaled 
time i it becomes exp(— p/^'"^) which converges to 
zero for T ^ at any fixed t > 1. 

We can now also come back to a point discussed at the 
end of Sec II V fjl in principle the function f{d) may not 
be strictly monotonically increasing, but could instead 
be constant from a certain value onwards, f{d) = foo 
for d > di,. Looking at (|45|l . however, this would lead 
to the conclusion that for any t < frx, the relaxation 
functions are nonzero, approaching a nonzero limit as 
i foo but then dropping discontinuously to zero (in the 
limit T ^ 0). This makes this hypothetical behaviour of 
f{d) rather unlikely. 

As explained, one of the main benefits of the super- 
domain model is that it allows one to access very long 
timescales that diverge extremely quickly as T decreases. 
Since this is precisely the regime that is hard to probe 
with simulations, comparing the theoretical predictions 
with numerical data is not straightforward. We chose 
e 0.02, the lowest value for which we can simulate a 
significant part of the equilibrium relaxation. As Fig. ^ 
shows, the relation between rescaled FPTs f and rescaled 



distances d = ed has not yet reached its e ^ form, so 
it would make no sense to compare numerical data with 
superdomain predictions based on the latter. Instead, we 
use the numerically obtained FPTs to link timescales and 
lengthscales. In detail, we obtain for each unnormalized 
distance d the measured FPT, and plot against this time 
the superdomain model's persistence function Vsdid) at 
the relevant scaled distance d — ed. Since we measured 
FPT in four different ways — for up- and down- flips, and 
averaging FPTs or log-FPTs — which we expect to be- 
come identical only for T — > 0, this procedure gives four 
slightly different curves for the superdomain predictions. 
As shown in Fig. |S1 these bracket the numerically simu- 
lated up-spin persistence and correlation functions rather 
well, especially given that there are no fit parameters in 
the comparison. The largest deviations occur for large 
times, corresponding to large lengthscales. This is con- 
sistent with Fig. 0] which shows that in this regime the 
behaviour still differs substantially from that expected in 
the hmit T ^ 0. 



V. CONCLUSION 

We have studied the dynamics of the East model. Con- 
sisting of uncoupled spins in a downward-pointing field, 
the model has trivial equilibrium statistics. However, 
the kinetic constraint that spins can only flip if their left 
neighbour is up causes pronounced glassy features in the 
dynamics at low T, when the concentration of up-spins 
is low. 

We first studied the non-equilibrium coarsening dy- 
namics after a quench to low T. In the limit e = 
exp(— 1/T) — > the equilibrium concentration of up- 
spins at the new temperature is negligible (« e) and the 
flipping-down of up-spins becomes irreversible to lead- 
ing order. This allows the dynamics to be described as 
coarsening via coalescence of down-spin domains. The 
process is hierarchical, being governed by a series of well- 
separated timescales. We solved this hierarchical coars- 
ening dynamics exactly, using an independent intervals 
method that becomes exact for T — > 0. Anomalous 
coarsening results, with typical domain lengths scaling 
as d ~ t"^'"^. The dominant divergence of the equi- 
libration time for low T can also be estimated, and is 
given by the factor = exp(l/r^ln2), an EITS depen- 
dence typical of fragile glasses. For large domain sizes d 
that are still small compared to the equilibrium value Jeq; 
the domain size distribution approaches a scaling form. 
We showed that this scaling distribution is equal to that 
of the paste-all model, and were able to define a whole 
family of interpolating models that all share this scaling 
distribution. 

In the second part of the paper we focussed on the 
equilibrium dynamics at low T. We showed that the 
standard relaxation functions, spin-autocorrelation and 
persistence of up- and down-spins, become identical for 
T ^ 0, so that only one of them needs to be consid- 
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ered. We then investigated the relation between time- 
and lengthscales. GeneraUzing from the results in the 
coarsening regime d = d/dcq ^ 1, we introduced a time 
scaling hypothesis. This implied that for T — > one 
has continuous timescale separation, with domains of any 
two different scaled sizes d relaxing on well-separated 
timescales. On this basis we proposed a model of super- 
domains, which are bounded by up-spins that are frozen 
on long timescales. The dynamics of this model is non- 
trivial and not, as yet, analytically tractable, but can 
easily be simulated since time is effectively measured in 
terms of the scaled distance d, whose relevant values re- 
main 0(1) even for T — > 0. We verified that the model 
obeys important consistency requirements, in particular 
the equality of up-spin and down-spin persistence, the 
fluctuation-dissipation theorem, and the stationarity of 
up-spin concentration and domain size distribution in 
equilibrium. 

From the superdomain model predictions we finally de- 
duced that the equilibrium relaxation functions should 
decay for low T as g{i), with g{-) a scaling function of 
i = (t/t*)-'"'"^. This demonstrates strong stretching for 
low T, but the overall relaxation is more complicated 
than a stretched exponential. In fact, the function g(-) 
decays faster than exponential, and in the limit T ^ Q 
at fixed t reaches zero at a finite value of t. The low- 
est temperature T that we can conveniently simulate, 
corresponding to e ^ 0.02, is still rather far from the 
asymptotic T — > limit but nevertheless showed rea- 
sonable agreement between numerical simulations and 



appropriately extracted predictions of the superdomain 
model. We would suggest that stretching but not simple 
stretched exponential behaviour may be rather generic 
in glassy dynamics. Actual stretched exponentials could 
more often than not be just convenient fitting functions 
over a limited number of decades in time. To clarify 
this point, a study of low-temperature relaxation in other 
solvable models exhibiting glassy dynamics would obvi- 
ously be desirable. 

In future work, it would be interesting to see whether 
the out-of-equilibrium response of spins to a local field 
could also be analysed within the irreversible coarsen- 
ing framework we used above. This response function 
was simulated in |^ and found there to be monotonic; 
at lower T, however, non-monotonicities should appear 
according to a later conjecture ^3d\. A closer investiga- 
tion of out-of-equilibrium FDT relations would also be 
worthwhile. Previous results for these have to be 
regarded with some caution since they were constructed 
using a disconnected correlator; see e.g. [s^ for discus- 
sion of this point. Finally, as regards the equilibrium 
dynamics, it will be interesting to analyse the implica- 
tions of the superdomain description for dynamic het- 
erogeneities, making a connection to the recent work of 
Garrahan and Chandler 21]. 
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